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Angular distributions of differential cross sections from the latest CLAS data sets I2J for 
the reaction 7 + p^K~^ + A have been analyzed using associated Legendre polynomials. 
This analysis is based upon theoretical calculations in Ref. ^ where all sixteen observables 
in kaon photoproduction can be classified into four Legendre classes. Each observable 
can be described by an expansion of associated Legendre polynomial functions. One of 
the questions to be addressed is how many associated Legendre polynomials are required 
to describe the data. In this preliminary analysis, we used data models with different 
numbers of associated Legendre polynomials. We then compared these models by calcu- 
lating posterior probabilities of the models. We found that the CLAS data set needs no 
more than four associated Legendre polynomials to describe the differential cross section 
data. In addition, we also show the extracted coefficients of the best model. 
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1. Introduction 

Significant information on the structure of the nucleon can be obtained by study- 
ing its excitation spectrum. Over the last few decades, a large amount information 
about the spectrum of the nucleon has been collected. Most of this information 
has been extracted from pion-induced and pion photoproduction reactions. How- 
ever, pionic reactions may have biased the information on the existence of certain 
resonances. Constituent quark model calculations predict a much richer resonance 
spectrum than has been observed in pion production experiments . Predicted res- 
onances which have not been observed are called "missing" resonances. Instead, the 
constituent quark model also predicts that these "missing" resonances may couple 
strongly to KA and KS channels or other final states involving vector mesons I SHE1 
Since performing kaon-hyperon, kaon-nucleon or hyperon-nucleon scattering exper- 
iments is a daunting task, kaon photoproduction on the nucleon appears to be a 
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good alternative solution 

Experiments on kaon photoproduction and electroproduction started in the 
1960s. However the old experimental data are often inconsistent and have large 
error bars. In recent years a large amount of data for kaon photoproduction has 
been collected. High statistics data from CLAS, for differential cross sections, recoil 
polarization, C x an d C* double polarizations for the reaction 7 + p — > K + + A 
have bee n publi shed Add itional experimental data have also been measured by 
saphirEEEoI LEPS lU112| and GRAAL [I3] 

Several previous a nalyses have been applied to the results of these ex periments, 
such as Isobar models |4 | 5 | 14|15 _U6| and Coupled channel models HZEHni However 
different theoretical model calculations often produce very different predictions. 

In Ref.l all sixteen observables in kaon photoproduction were shown to be classi- 
fied into the classes £o(I; E; C z /; W), £i a (P; H; C x /; L x /), £15 (T; F; O x /; T z /) and 
£2(^5 G; O z /; T x /), where each class is an expansion in a different set of associated 
Legendre polynomials. What is not apparent is how many terms in each expansion 
are required. This work attempts to address the issue by examining data models 
with different numbers of terms, and calculating which one has the greatest posterior 
probability. In this article we only focus on the differential cross section observables, 
which are described by the associated Legendre class £o- 

2. Analysis Procedure 
2.1. Data Model 

We construct data models based on Legendre class £q. These data models can be 
written compactly as follows: 

l=L 

M l° =Y. AlPl ° ( cos °)' W 

1=0 

where m£° is the data model, and A\ and Pw(cosO) are the coefficients and asso- 
ciated Legendre polynomials. Each data model therefore has a different "order" or 
maximum number of polynomials. Our task is to find the most likely order. 



2.2. Model Comparison 

To determine the best model, we evaluate the posterior probability ^Q] f or eac ] 1 
data model. The ratio of the probabilities for Ml and Mo can be written, using 
Bayes theorem, as follows : 

_ P(M L \D) = P(D\M L ) P(M L ) 

P(M \D) P(D\M ) P(M ) ' 1 ] 

where P(Ml\D) is the posterior for the Ml model, P(D\Ml) is the probability 
that the data would be obtained, assuming Ml to be true (the likelihood). With 
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no prior prejudice as to which variant is correct, we obtain the ratio of likelihoods: 

_ P(D\M L ) 
P(D\M ) 

The likelihood P(D\M L ) is an integral over the joint likelihood P(D,{Ai}\M L ), 
where {Ai} represents a set of free parameters: 

P(D\M L ) = J ... J ' P{D,{Ai}\M L )d L Ai, 

= j ... j P(D\{Ai},M L )P({Ai}\M L )d L Ai. (3) 

The function P(Ai\Ml) is the prior probability that the parameters take on specific 
values. We assume that each parameter A\ lies in the range A™ %n < A\ < A™ ax ', and 
we can write the prior as the reciprocal of the volume of a hypercube in parameter 
search space as P({Ai}\Ml) = L ] x — —. — . If the errors in the data points 

11 K A i ~ A i ) 

are Gaussian, it can be shown that P(D\{Ai} 1 Ml) oc exp ^— where x 2 is the 
sum of squared residuals. Using a Taylor series expansion about the minimum x 2 ? 
X 2 ~ Xmin + |(-^- — X ) T V 2 x 2 (X — X ) + we can write an approximate form 
for the likelihood: 

P(D\M L ) oc — T L\(^) L fxLA (4) 

[] L (A™ ax - A™ n ) x ^Det(VV X 2 ) \ 2 / 

where L is the dimension of the integral and (Det(\7Vx 2 ) is the determinant of the 
Hessian matrix, which in turn is the inverse of the covariance matrix. 



3. Results 

Using the above analysis procedure, for each of the available photon energy bins, we 
fitted each data model to the angular distribution. This was carried out using the 
standard minimization package MINUIT. We then compared models with different 
numbers of Legendre polynomials by evaluating Eq. (|4]) for each data model. 

To illustrate the procedure, we first choose one photon energy bin at E 1 = 1.824 
GeV as an example. The posterior probabilities are shown in Fig.[lJ where the order 
of the data model is shown on the horizontal axis. The maximum posterior is given 
by the data model containing four associated Legendre polynomials. 

On the left side in Fig. [2] we show the fit of the fourth order data model to 
the CLAS differential cross section data ^ for E 1 = 1.824 GeV. The procedure 
is repeated for each photon energy bin. The right side in Fig. [2] the order of data 
model which has the greatest probability at each photon energy is plotted. It can be 
seen that this generally increase from threshold into the resonance region, but that 
the maximum is mostly at the fourth order. The distributions of the polynomial 
coefficients for fourth order data models as a function of photon energy is shown in 
Fig. El 
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Fig. 1. The posterior probabilities for different orders of data model, for E 1 = 1.824 GeV. 

4. Conclusion 

We have analyzed the Legendre polynomial decomposition of differential cross sec- 
tion data. We generated data models with different numbers of associated Legendre 
polynomials. We then compared them by calculating posterior probabilities. From 
this analysis, we found that differential cross section data in this case requires at 
least four associated Legendre polynomials. 
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Fig. 2. Plot showing the fit of the fourth order data model to the CLAS cross section data for 
E 1 = 1.824 GeV (on the top). The order of associated Legendre polynomial for all photon energy 

(on the bottom). 
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Fig. 3. Extracted associated Legendre polynomial coefficients for each photon energy. 



